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Abstract 

Building on the previous results of the Weyl chamber steering method, we demonstrate how to 
generate high-fidelity controlled-not by direct application of certain, physically relevant Hamiltonians 
with fixed coupling constants containing Rabi terms. Such Hamiltonians are often used to describe 
two superconducting qubits driven by local r/-pulses. 

It is found that in order to achieve 100% fidelity in a system with capacitive coupling of strength 
g one Rabi term suffices. We give the exact values of the physical parameters needed to implement 
such CNOT. The gate time and all possible Rabi frequencies are found to be t = n/(2g) and Sl\/g = 
V64n 2 — 1, n = 1,2,3,.... Generation of a perfect CNOT in a system with inductive coupling, 
characterized by additional constant k, requires the presence of both Rabi terms. The gate time is 
again t = ir/(2g), but now there is an infinite number of solutions, each of which is valid in a certain 
range of k and is characterized by a pair of positive integers (n, m), 



We distinguish two cases, depending on the sign of the coupling constant: (1) the antiferromagnetic 

case (k > 0) with n > m = 0, 1, 2, . . .; (2) the ferromagnetic case (k < 0) with n > m = 0, 1, 2, 

The paper concludes with consideration of fidelity degradation by switching to resonance. Simu- 
lation of time-evolution based on the 4th order Magnus expansion reveals characteristics of the gate 
similar to those found in the exact case, with slightly shorter gate time and shifted values of Rabi 
frequencies. 

PACS number (s): 03.67.Lx, 03.65.Fd, 85.25.-j 
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1 Introduction 



The goal of this paper is to show how to successfully generate, with a minimal amount of effort, a 
high-fidelity CNOT in systems consisting of coupled superconducting qubits described by Hamiltonians 
with coupling constants fixed by architecture. We want to find out if making a CNOT in such systems 
can, at least theoretically, be reduced to a simple push of a button. Once initially set, the system 
itself would then take care of the needed gate by freely evolving from the identity of the unitary group 
to the target along the corresponding geodesic. 

We thus approach the subject of time-control from the "opposite", non-optimal end. The plug- 
and-play implementations considered in this paper may be useful for simple gate designs which do 
not rely on sophisticated time-control of their physical parameters. The question of simplicity will 
also become important once the best qubit is identified, and put together with thousands of others to 
form an integrated circuit. Optimally controlling [lj qubits in such a circuit would become a difficult 
task, potentially leading to a compromise between the optimal and the simple. Thus, in this paper we 
limit our attention mostly to time-independent Hamiltonians. The only place where time dependence 
necessarily arises is in the analysis of fidelity degradation due to switching to resonance, considered 
in Section El 

Two important developments motivated our work. The first was a series of experiments [21 [3j H] 
which showed that the macroscopic quantum states [51 E] of Josephson junctions have coherence times 
long enough for such states to be used as qubits. That led to our choice of the physical Hamiltonians. 
The other development was the work on the Weyl chamber steering method [3 [8] in which several 
examples of direct generation of various important gates, including CNOT, was considered from a 
purely geometrical standpoint. 

In what follows we will use the square brackets [• • •] to denote the gate — a representative of a 
given local class — located in a Weyl chamber. 

2 The concept of local equivalence 

We begin by recalling the notion of local equivalence within the special unitary group ££7(4). Two 
gates U\,Vi £ ££7(4) are locally equivalent if they differ only by local transformations: 

U x ~ Vi ^ Ui = h Vi k 2 , ki,k2 € ££7(2) (g) ££7(2). (1) 

It is obvious that local equivalence is an example of equivalence relation: it is reflexive (£7 ~ £7), 
symmetric (£7 ~ V -4=>- V ~ £7), and transitive (£7 ~ V, V ~ W =^ £7 ~ W). When an equivalence 
relation is given on a set, it partitions the set into disjoint classes. Thus the entire group ££7(4) is 
partitioned into classes of locally equivalent elements; any such element (a two qubit gate) belongs to 
one and only one local equivalence class. 

To generalize the concept of local equivalence to the full £7 (4) we use the fact that any £7 G £7(4) can 
be written in the form U = e ia U\, with U\ G ££7(4) (even though, as groups, £7(4) ^ £7(1) (g> ££7(4)). 
Then two gates £7 = e ia U±, V = e 1 ^ V\ G £7(4) are locally equivalent if their representatives in 5£7(4) 
are equivalent: 

U ~V <^=>f/i~Fi. (2) 

One way to represent local classes in S'£7(4) is with the help of a Cartan decomposition of the 
underlying Lie algebra L = su(4). This can be seen as follows: 

For any semi-simple Lie algebra L, its Cartan decomposition is a direct sum of two subspaces 
L = K ® P, one of which, K, is itself a Lie algebra, [K, K] C K, and the other, P, satisfies 
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[P, P] C K, [P, K] C P. Because K n P = 0, any subalgebra A C P must necessarily be abelian. 
Such maximal abelian subalgebra Ac C P is called the Cartan subalgebra of L relative to the given 
Cartan decomposition. An important Cartan decomposition theorem (for groups) states that if a 
semisimple algebra L has decomposition L = K © P, then any element U± E exp L can always be 
written in the form 

Ui = h V 1 k 2 , (3) 

where V\ E exp Ac-, and E exp K. 

An example of a Cartan decomposition is provided by su(2), for which we can take K = span{|o" z }, 
P = span{!<r.j;, \o-y), and Ac = span{^a x } . 

Cartan decompositions of sit (4) are more complicated, but a special choice pQ of K spanned by 
the local operators 

x i = 2°i> Yl = 2 a y> Zl = 2°"*' ^ 2 = 5°"*' 12 = 2^' ^ 2 = 5°"*' ^ 

and P spanned by all nonlocal operators | o"qCj| (a, = x,y, z) is particularly useful. If we now take 
the span of mutually commuting operators 

XX= l -o- 1 x a 2 x , YY= l -ay y , ZZ = % -a\al (5) 

to serve as Ac, the Cartan decomposition theorem applied to SU(A) = exp[su(4)] gives 

Ui = k x exp (ciXX + c 2 yy + c 3 ZZ) k 2 . (6) 

This results in a convenient representation of SU(4:) local equivalence classes by the class vectors 
[ci, 02,03] E Ac- Representation of the U{4) classes is up-to-phase the same, 

U = e ia ki exp {c x XX + c 2 YY + c 3 ZZ) k 2 . (7) 

The correspondence between class vectors and equivalence classes is not one-to-one: there are 
certain symmetries which map class vectors to other class vectors of the same equivalence class. 
However it was shown in [7] that the correspondence can be made unique by restricting class vectors 
to a tetrahedral region in Ac, called a Weyl chamber. One such chamber, denoted by a + , is chosen 
to be canonical — it is described by the following three conditions [3 [9]: 

• 7T > C\ > c 2 > C3 > 0, 

• C\ + C 2 < 7T, 

• if C3 = 0, then c\ < ir/2. 

Even though the choice of the Weyl chamber is not unique, when fixed, it becomes a powerful tool for 
analyzing time evolution on the full U(4). 

In the original papers on the subject [8], several implementations of various useful gates by 
steering on a + have been demonstrated. Of particular interest to us is the way [CNOT] = [n/2, 0, 0] E 
a + was generated by a single application of Hamiltonians with controllable Ising type interaction, such 
as 

H vy = flixo-l + VL 2x a 2 x + VLi z o l z + VL 2z a 2 z + go\a\, (8) 

and similarly for H xx , H zz . In [3 [8], the local invariants of the corresponding unitary evolution U(t) = 
ex.p(iH yy t) were first found in closed analytic form. Numerical analysis then produced the values for 
the time and the coupling constants needed to achieve a perfect [CNOT] in just one application. The 
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yy and zz cases are particularly interesting because their interaction parts (when viewed as elements 
of Aq) point in the direction "perpendicular" to [CNOT]. 

Let us now turn to Hamiltonians used to describe coupled superconducting qubits. Such Hamilto- 
nians may contain several Ising terms, which we choose to be fixed by architecture (c/. \10\ II 1| ) . Our 
goal is to find out if they are also capable of generating [CNOT] in a single application, and if so, how 
many Rabi terms are required in order for such generation to be successful. 

3 Hamiltonians for coupled superconducting qubits 

When restricted to the two-qubit subspace, the Hamiltonian for two capacitively coupled phase qubits 
driven by rf-pulses is given by 

2 2 

tffast = -^^ + ^acos(c4t + <5iK + 5 ^, g>0. (9) 

i=l 1 i=l 



H V 

The level spacings are tunable between about 7 to 10 GHz. The Rabi frequencies 0^ and the phases 
8i are fully controllable. 

As it stands, this Hamiltonian is not suitable for actual gate implementations because of the very 
small time scale (t sc ~ 0.1 ns) associated with a^. The scale roughly corresponds to the time it takes 
to make a "round trip" on SU(4:) when going in the direction of — %H G su(4). To attain high fidelity 
of the gate we would have to control the pulses with accuracy which is experimentally unattainable. 

To avoid this problem, the rotating wave approximation (RWA) is used. We set the frequency of 
the external pulse to be equal to the transition frequency of the system (resonance), = u>i = u, 
bring the phases to zero, and work in the interaction picture, in which the time evolution is governed 
by 



H 



iH ty e -iH t 
2 



- sin(2cL , t') 
Oj cos 2 {ujt)a l x H a y +g cos(ujt)<jy — sm(ivt)al cos(wt)<T^ — sm(Lut)a 2 



i=l 



(10) 



After averaging over fast oscillations the Hamiltonian becomes 



or in the language of Lie algebra su(A), 

- %Ht = -t [J2i Xi + Q 2 X 2 + g (XX + YY)} . (12) 

There is an interesting modification of this qubit architecture based on inductive coupling between 
the qubits (see Appendix). In that case an additional coupling constant k is introduced resulting in 
the full Hamiltonian 

- %Ht = -t [Oi X x + tt 2 X 2 + g {XX + YY + kZZ)] . (13) 
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4 Implementing [CNOT] by time-independent Hamilto- 
nians 



4.1 The XX + YY case 

Here we show how to implement the perfect [CNOT] by a single application of the experimentally 
available Hamiltonian (1121) with one Rabi term 



iHt = -t [ft], Xi + g (XX + YY)] 



(14) 



We first establish a useful group-theoretical result about the exponentials of (|14p . We consider a 
subalgebra of sw(4), ^4o = span{Xi, XX, FY, ZY}, whose generators obey the commutation relations 
summarized in the following table: 



Ac 





Xi 


XX 


YY 


ZY 


Xi 








-ZY 


YY 


XX 














YY 


ZY 








-x x 


ZY 


-YY 





Xi 






(15) 



Notice that ZY generator had to be added to the set of generators in order to ensure closure under 
commutation. 

The generator XX is central in Aq. Its span is a one-dimensional abelian ideal I C Aq, which makes 
Aq into a nem-semisimple Lie algebra. This algebra is a direct sum of two subspaces, Aq = I© Ai, 
with A\ = span{Xi, YY, ZY} isomorphic to su(2). 

Exponentiation of Aq results in a /oitr-dimensional subgroup Go = exp^o C SU(4) whose every 
element W can be written product 



W _ e -ci XX e -a X! -c 2 YY —b ZY _ 
an element of exp A\ 



(16) 



Because A\ has its own Cartan decomposition (e.g., A\ = span{Xi} © sp&n{YY, ZY} with Cartan 
subalgebra Ac = span{YY}), application of the Cartan decomposition theorem to the group G\ = 
exp A\ gives 

W = e -ciXX e -aX le -c 2 YY e -pX K 

Now, because e~ lHt = e-t^x^+gtyXX+YY)} [ s itself an element of Go = expAo, it can also be 
written as such product 



e -t[n 1 X 1 +g(XX+YY)] _ g-aXi g - Cl XX-c 2 YY e ~/3X 1 



(18) 



where the commutativity of X\ and XX has been used to rearrange the factors on the right. 

Thus given any gate (t,£l\,g) generated by H in (I14j) . with coupling g fixed by architecture, there 
always exists a set of four numbers (a, c\, C2, ft) which represent that gate exactly. This is essentially 
the same good old Cartan decomposition theorem ([6]) for SU(4), restricted to a physically relevant 
subgroup Go associated with experimentally available Hamiltonian (JHJ). 

By fixing J7i and evolving the system freely under the action of its Hamiltonian, we generate a flow 
along a straight line in sit (4) which has its image in both, the Weyl chamber and the group SU(4). 

The converse is unfortunately not necessarily true: for a given gate (a, c\, C2, ft) it is not always 
possible to find (t, fii) generating that gate. This is because of noncommutativity between X\ and 
YY. 
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Even though there is no guarantee that H would be able to reach every (ci, C2, 0) equivalence class 
in a single application, it is interesting to find out what actually can be reached. In particular, we 
want to know if [CNOT] is reachable. 

It is obvious that if Q,\ = 0, then any local class with c\ = C2 and C3 = can be directly reached. 
What if Oi ^ 0? 

We will be interested in purely geodesic solutions to (fl~8j) . Thus, setting a = (3 = and directly 
exponentiating in (I18p we get the matrix equation 



cos b± — cos a 


icd sin a 
z(sin b\ — d sin a) 





cos 62 — cos a 
i (sin 62 + dsina) 
icd sin a 



icdsina i(sin 61 — <isin a) 

i (sin 62 + (ism a) icd sin a 

cos 62 — cos a 

cos 61 



cos a 



0, 



where 




'1,2 



-(ig - ci ± c 2 ) 



If f^i 7^ 0, then sin a = 0, leading to cos a 
possible if b\ = 62 + 27rm, or 

C2 = 2irm, 
The minimal time solution is therefore 

t ■ = - 1 
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±1 and cos 61 = 
0,1,2, 



cos 62 



±1, sin 61,2 



(19) 



(20) 



m 



0. This is only 
(21) 




n 



1,2,3,.... 



(22) 



Thus, by following H along its geodesic for a time t m ; n we can (exactly) reach any equivalence class 
located on the XX axis of the Weyl chamber. It is clear that the coupling constant g plays the role 
of the scaling factor between the coordinate c\ of the point on this axis and the actual, physical time 
required to reach that point. 

We want to emphasize that the corresponding trajectory in the Weyl chamber is not a straight 
line from O to [CNOT], even though is is a straight line in the full Lie algebra su(4). This is because 
with fti now being fixed by (J22J), it is not true that e -t[^i+9(xx+YY)] ^ e - Cl (t)xx ig sat i s fi ec i f or 
all intermediate times < t < i m in- In general for such times, C2(t) / in (fl8|) . 

One important local class belonging to the XX axis is [CNOT], which can be made by expo- 
nentiating the Hamiltonian — iH^Q^Q^t = —tXX with t = tt/2. The corresponding matrix in the 
computational basis is 



[CNOT] = exp (--XX 



7T 





1 








— i 


1 





1 


—i 





72 





—% 


1 







—i 








1 



6 517(4). 



(23) 



The Makhlin invariants [12] of this [CNOT] are G\ = and G2 = 1, as required, so this is a perfectly 
acceptable [CNOT] out of which the canonical CNOT with matrix 



CNOT 



eU(A) 



(24) 
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can be produced by local manipulation of individual qubits. We will give some examples of the needed 
for this local r/-pulses in Section [5j 

Thus, setting c\ = tt/2 in ([22]) . we get 

imin = 7r/(2s0, 



Qi/g = V64n 2 - 1, n = l,2,3,.... (25) 
The simplest exact [CNOT] is then 

[CNOT] = e"t IVm 9 x 1+9 (xx + yy)] (2g) 

For completeness we mention another [CNOT] solution to (|18|) . 



U x /g = Vl6n 2 - 1, n = 1,2,3,... , (27) 

which is up-to-phase geodesic, 

[CNOT] = e M e-£ ^x 1+g (xx+YY)] (28) 



4.2 The XX + YY + kZZ case 

We are now going to show that a single application of the Hamiltonian given in (|13p generates the 
exact [CNOT], provided both Rabi terms are kept. 

Derivation parallels the previous case. We concentrate on subalgebra Aq C su(4): 



Aq 





Xx 


x 2 


XX 


YY 


ZZ 


YZ 


ZY 


x x 











-ZY 


YZ 


-ZZ 


YY 


x 2 











-YZ 


ZY 


YY 


-ZZ 


XX 























YY 


ZY 


YZ 











-X 2 


-Xx 


ZZ 


-YZ 


-ZY 











Xx 


x 2 


YZ 


ZZ 


-YY 





X x 


-X x 








ZY 


-YY 


ZZ 





Xt 


-x 2 









(29) 



The generator XX is again central in Aq, spanning a one-dimensional abelian ideal / C Aq. This 
makes Aq into a non-semisimple Lie algebra which is a direct sum of two subspaces, Aq = I © Ax, 
where A 1 = span-pfi, X 2 , YY, ZZ, YZ, ZY}. 

Exponentiation of Aq results in a seven-dimensional subgroup Go = exp^o C S'C/(4) whose every 
element W can be written as a product 



w _ e ~ciXX g— 0X1—6X2— C2 YY—C3 ZZ—dYZ—fZY 
an clement of exp A\ 



(30) 



Because A\ has its own Cartan decomposition (e.g., Ax = span{Xi,X2} © span{YY, ZZ,YZ, ZY} 
with a Cartan subalgebra Ac = spanjyY, ZZ}), application of the Cartan decomposition theorem to 
the group Gx = exp A\ gives 



w = e -ciXX e - a X 1 -pX 2 e -c 2 YY-c 3 ZZ e - v X 1 -^X 2 ^ 



(31) 
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Since e~ iHt = e -miXi+n 2 x 2 -g(xx+YY+kZZ)] ig itgelf an e i ement f q q = exp A , it can also be 
written as such product 

e -t[Q. 1 X 1 +Q. 2 X 2 +g{XX+YY+kZZ)} _ g-aXj-^Xj g-ci XX-c 2 YY-c 3 ZZ ^X^Xi ^) 

Thus, given any gate (t, fix, $1%, g, k) generated by H in (fl~3l) . with couplings g and k fixed by 
architecture, there always exists a set of seven numbers (a, /3, ci, C2, C3, 77, £) which represent that 
gate exactly. This is again the Cartan decomposition theorem for SU(4), this time restricted to a 
physically relevant subgroup Go associated with experimentally available Hamiltonian (|13p . 

With a = /3 = r/ = ^ = 0, equation (|32[) has an infinite number of minimal time solutions labelled 
by pairs of integers (n, m) and valid in certain intervals of k, 

t ■ = ^ 

g 



o, 2 _ (^y.^p (^y-(t±i)\ 
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For [CNOT], (ci,c 2 ,c 3 ) = (vr/2,0,0), which gives 



l^-( h -^) ±Jie^-( h -±l) . (34) 



, v-- v 2 / ± V 16m2 -B 1 " 

Two cases will be distinguished here, depending on the sign of the coupling constant. 



4.2.1 kZZ coupling of antiferromagnetic type 

In this case, k > and n > m = 0, 1, 2, 3, The previous result for = is correctly recovered when 

n = m = 0. For n = m = 1, 2, . . . we get an (n, n)-th solution valid in the interval < k < 8n — 1. 

Below, we list Rabi frequencies required to generate a perfect plug-and-play [CNOT] in the anti- 
ferromagnetic case for several values of n and m: 



k 


n^/g n^/g 


^' 1} /g n^/g 


^/g ^/g 


^ S) /g 


0.00 
0.10 
0.25 
0.50 
1.00 
3.00 
5.00 
7.00 


7.937254 0.000000 
7.936614 0.012600 
7.933253 0.031513 
7.921238 0.063121 
7.872983 0.127017 
7.337085 0.408882 
6.109853 0.818350 
2.645751 2.645751 


11.952987 4.015733 
11.949341 4.025327 
11.942076 4.040336 
11.925151 4.067034 
11.872983 4.127017 
11.401356 4.473152 
10.391718 5.100215 
7.416198 7.416198 


15.968719 0.000000 
15.968405 0.006262 
15.966755 0.015658 
15.960859 0.031327 
15.937254 0.062746 
15.683221 0.191287 
15.162165 0.329768 
14.344402 0.487995 


23.979157 0.000000 
23.978949 0.004170 
23.977852 0.010426 
23.973935 0.020856 
23.958261 0.041739 
23.790420 0.126101 
23.451110 0.213209 
22.932659 0.305242 


9.00 
11.00 
13.00 
15.00 






13.173201 0.683205 
11.536501 0.953495 
9.164486 1.418519 
3.872983 3.872983 


22.222421 0.404996 
21.301017 0.516407 
20.139099 0.645511 
18.691066 0.802522 


17.00 
19.00 
21.00 
23.00 








16.881526 1.007018 
14.570504 1.304004 
11.429081 1.837418 
4.795832 4.795832 



(35) 



S 



Notice that the product of Rabi frequencies for solutions of (n, ra)-type is always equal to k, 

n (n,n) n (n,n) /g2 =k _ (36) 

4.2.2 kZZ coupling of ferromagnetic type 

In this case, k < and n > m = 0, 1, 2, A particularly interesting family of solutions occurs for 

m = 0. These solutions (c/. ([27]) ) exist only for the coupling constant k = — 1: 

fi^ 0) /5 = Vl6n 2 - 1, ra= 1,2,3,.... (37) 
For certain n > m = 1, 2, 3 we find: 



k 


n^/ 9 




n 


S"»/ 9 


n^/ 9 






n nn 

— u.uu 


11.952987 


4 


015733 


15 


958206 


8.020952 


1 n ovqoqo 

ly.y i ovot) 


4.uuoziy 


-0.10 


11.955678 


4 


006464 


15 


961996 


8.012782 


19.974723 


4.000055 


-0.25 


11.957932 


3 


993165 


15 


966096 


8.001330 


19.974919 


3.992507 


-0.50 


11.956946 


3 


972586 


15 


968719 


7.984360 


19.972632 


3.980447 


-1.00 


11.937254 


3 


937254 


15 


958261 


7.958261 


19.958261 


3.958261 


-2.00 


11.826744 


3 


889490 


15 


874508 


7.937254 


19.890241 


3.921521 


-3.00 


11.618950 


3 


872983 


15 


705143 


7.959176 


19.769413 


3.894906 


-4.00 


11.307441 


3 


891243 


15 


444795 


8.028595 


19.594811 


3.878578 


-5.00 


10.880300 


3 


952097 


15 


083052 


8.154848 


19.364917 


3.872983 


-6.00 


10.316246 


4 


071248 


14 


600739 


8.355741 


19.077582 


3.878898 


-7.00 


9.573955 


4 


282452 


13 


959460 


8.667957 


18.729907 


3.897510 


-8.00 


8.550870 


4 


677887 


13 


060789 


9.187806 


18.318045 


3.930550 


-9.00 


6.244998 


6 


244998 


10 


908712 


10.908712 


17.836915 


3.980509 


-11.00 
-13.00 
-15.00 
-17.00 






16.637303 
15.038297 
12.817255 
7.937254 


4.147307 
4.455292 
5.071289 
7.937254 



(38) 



5 Local r/-pulses 

Here we give two (ideal) local r/-pulse sequences needed to produce the true, canonical CNOT from 
the [CNOT] located in the Weyl chamber via 

CNOT = i 77/4 ki [CNOT] k 2 . (39) 

Sequence 1: Total rotation ip = 3tt/2, 

fel = e -fn e |(x 1 -x 2 ) ) fc 2 = e fn (40) 

Sequence 2: Total rotation (p = 3tt/2, 

k x = e-^ Y \ k 2 = e~^ Zl e^ {Xl ~ x ' 2 \ (41) 
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6 Fidelity degradation by switching to resonance 

We are now going to consider the effect of switching — the actual process of bringing the qubits to 
resonance with external r/- field. 

We will assume that [CNOT] generation is performed in three stages: 

Stage 1. Tuning: Very fast linear rise of Rabi frequencies to their optimal values (for a given k), 
with time-dependent Hamiltonian 



tunc 



(n° pt Xi + n° pt x 2 ) —!— + g{xx + yy + kzz) 



(42) 



where e <C 1 characterizes the fraction of the total gate time spent in the switching mode 
(tune/detune). We denote the corresponding time evolution by ?7 tune ; 

Stage 2. Resonance: Evolution with optimal Rabi frequencies for most part of [CNOT] generation, 

^resonance; where 



-iff, 



resonance 



resonance 



ft° pt Xi + n 2 pt x 2 + g{xx + YY + kZZ) 



(43) 



and 



Stage 3. Detuning: Fast linear decrease of Rabi frequencies to zero, with time-dependent Hamil- 
tonian 



%H, 



detune 



(n opt x 1 + n° pt x 2 ) [i 



t 



et 



+ g{XX + YY + kZZ) 



:ate 



The full gate is then the product 

[CNOT] t — ^detune ^resonance ^tune • 



(44) 



(45) 



In order to find the optimal values of the gate parameters, we first truncate the Magnus expansion 
for each switching part of the gate: 



^switching (f^gate) — 6 



c(e igato) 



with 



a (etgate) 



gate 

~1> I ^ -^switching (^) 

10 

.^2 reta 



(46) 



(47) 



{-iy 



i gate 



dt 



-^switching (^1 ) i / dt 2 -^switching (^2 / 
JO 



and then use the 4th order iteration method [13J 

U n+1 = e an U n , 
h 



(48) 



Oft 

A 2 
h 



- (A, + A 2 ) + [A 1} A 2 ] 

i-ff switching (^n Cl/l), 

iH switching (tn C 2 h), 
£ ^gate 



N ' 
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where N is the number of iteration steps and 



Cl,2 



1 

2 3 



(49) 



are the nodes of the Gauss-Legandre 4th order quadrature in [0, 1]. As described in [15] ; despite the 
presence of only /i and /i 2 terms, this is indeed the 4th order approximation to the Magnus expansion. 

Optimization of gate's parameters was performed using Nelder-Mead simplex direct search with 
bound constraints p3] for the minimum of the square of the Frobenius distance 



d 



Frobenius 



tr 



[CNOT] - [CNOT] opt ) ([CNOT] - [CNOT] 



opt 



between the optimal [CNOT] opt given by (|45p and the perfect [CNOT] given in 
enormous richness of possible CNOTs, we optimized only the most interesting of them. With N 
we simulated the switching mechanism for e = 0.025: 



(50) 

Due to the 
2500 
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*gate, X7r/(2s0 


<7 2 v i n~ 3 

"Frobenius J A lu 


0.50 


8.130446 


0.064667 


0.999268 


1.9 


0.25 


8.141971 


0.032287 


0.999363 


1.6 


0.10 


8.145193 


0.012910 


0.999390 


1.5 


0.00 


8.145807 


0.000000 


0.999395 


1.5 


-0.10 


12.269887 


4.111716 


0.999382 


1.5 


-0.25 


12.272602 


4.098223 


0.999348 


1.6 


-0.50 


12.272970 


4.077597 


0.999237 


1.9 



(51) 



Using the local pulses given in Section [5] we can transform the corresponding optimal [CNOT] opt 
to its canonical form. For example, 



CNOT fe=0 .25 



0.9998 - 0.0025i -0.0001 - 0.0017f 

-0.0001 - 0.0017?' 0.9998 - 0.0025? 

0.0193?' 0.0047?' 
-0.0001 - 0.0047? -0.0192? 



-0.0001 - 0.0047? 
-0.0193? 



0.0192?' 
0.0001 + 0.0047?' 



0.0001 - 0.0007?' 0.9998 + 0.0025?' 
0.9998 + 0.0025?' 0.0001 - 0.0007?' 



whose intrinsic fidelity relative to the computational basis is 



1 4 

F = -J2\ (CNOT 1 • CNOT fc=0 . 



2". 



r=l 



0.9998. 



(52) 
(53) 



7 Conclusion 

To recapitulate, we have demonstrated how to implement various high-fidelity, plug-and-play CNOT 
logic gates by a single application of Hamiltonians used to describe coupled superconducting qubits. 
Architectures involving capacitive and inductive couplings have been analyzed for which the physical 
parameters required to generate the perfect CNOT have been found in closed, analytic form. Ad- 
ditional numerical simulations based on 4th order Magnus expansion were used to correct this ideal 
limit by taking into account the switching to resonance. 
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Appendix: Derivation of the Hamiltonian for inductively 
coupled flux qubits 

The system consists of two superconducting loops of self-inductance L coupled by a mutual inductance 
M and driven by external magnetic fluxes <3?j, as shown on Fig. 1. Each loop is interrupted by a 
Josephson junction of critical current Iq and effective capacitance C. 

Inside the superconducting material, the Ginzburg-Landau complex order parameter can be written 
in the usual form, 

^)e ie{r \ (A-l) 



^(r) 

where p s is the density of Cooper pairs. The supercurrent flux is 

m* V n 1 



(A-2) 



where A s is the vector potential inside the superconductor and q = — 2e. The gauge is assumed to be 
chosen in such a way as to guarantee n • j s = at the surface. In the clean limit, j s = inside the 
material, which leads to 

A s = ~V0. (A-3) 

q 

The line integral of A around Loop 1 is (here, the limits of integration 1 and 2 denote the two 
sides of the junction, the rest of notation is clear from context): 
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-e x )- f A s d\ 
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r r 2 










01 - 2 + V0cfl 






q 







n + f V0dl 

'Loop 1 



!>i + 2irn 



(A-4) 



where (pi is the gauge invariant phase difference across the first junction, and n\ = 0, 1, 2, 3, . . .. This 
"quantization" condition is justified by the fact that the order parameter is nonzero inside the junction 
J. Alternatively, 



Loop 1 



Adl = $| xternal + $f lf + ^mutual = $1 _ LIl + M J 2 _ 



(A-5) 
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Introducing 

ft <£ sc ^ 
a = — = —, Ej = al , uj 

2e Z7T 



'LC 



a 2 hW (2e) 2 
Eq = — = , &c ~ 



L 2E C 



2C 



T = f, (A-6) 



and taking into account (|A-4|) and (|A-5|) together with the Josephson equations, we get the equations 
of motion for the two gauge-invariant phase differences, 



-Ej (sin 0i — T sin 02) — Eq I 0i + 2im 



a 2 C (0i - T0 2 

a 2 C (J>2 — T<pij = —Ej (sin 02 — T sin 0i) — -Eo ( 02 + 27m 2 



2tt$i 



(A-7) 



It is clear that 2-rrni are not dynamical and can be absorbed into 0j + 2irni — > 0j without affecting 
any physics. Notice also that in the limit M — > we correctly recover a well-known result for two 
independent loops. 

Because of the presence of terms T i?j sin 0j in (|A-7|) . we cannot immediately deduce system's 
Lagrangian. However, by multiplying the second equation in (lA-7h by T, and by adding the result to 
the first equation, we can separate the variables 



O' 



'■C (l - T 2 ) 0i = -Ej (l - T 2 ) 51110! - Eq 



2vr$i\ / 27T* 5 

n ^ + 1 02 -r — 



(A-6 



and similarly for 02. The Lagrangian is now easily found to be 

„2 



L 



(l - T 2 ) | ^ (0i 2 + 2 2 ) + Ej(co S 0i + cos 2 ) 



Eg 
2 



27T$i x 2 



+ 



2tt<I>.- x 2 



<Pl 



2vr$i 
"ST 



2tt$ ; 



5>s 



In the limit of weak coupling (T <C 1) the terms ~ T 2 can be dropped, giving 

i? = ^t + C/(0 1) 2 ), 



2m 



with = m0j, m = a 2 C, and 

2 



C/(0i,0 2 )=E 



i=l L 



E J ccek + T U i - — 



+ TE 



91 



27r$] 

$sr 



After projecting onto the qubit subspace and using RWA we recover (fT3 



2TT^r 



(A-9) 



(A-10) 



(A-ll) 
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